Description of the fluctuating colloid-polymer interface 
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To describe the full spectrum of surface fluctuations of the interface between phase-separated 
colloid-polymer mixtures from low scattering vector q (classical capillary wave theory) to high q 
(bulk-like fluctuations), one must take account of the interface's bending rigidity. We find that the 
bending rigidity is negative and that on approach to the critical point it vanishes proportionally to 
the interfacial tension. Both features are in agreement with Monte Carlo simulations. 
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One of the outstanding theoretical problems 
in the understanding of the structure of a simple 
liquid surface is the description of the full spec- 
trum of surface fluctuations obtained in light 

ulations [3|, y, |5|]. Insight into the structure of 
a simple liquid surface is provided by molecu- 
lar theories [0, LZ|], such as the van der Waals 
squared-gradient model, on the one hand and 

n Q 

the capillary wave model [8|, |9| on the other 
hand. The theoretical challenge is to incorpo- 
rate both theories and to describe the spectrum 
of fluctuations of a liquid surface from the molec- 
ular scale to the scale of capillary waves. 

Here, we report on a theoretical description 

n 

of Monte Carlo (MC) simulations 4] of a system 
consisting of a mixture of colloidal particles with 
diameter d and polymers with a radius of gyra- 
tion i?g. The presence of polymer induces a de- 
pletion attraction [10| between the colloidal par- 
ticles which may ultimately induce phase sepa- 
ration [ll|, [l^. The resulting interface of the 
demixed colloid-polymer system is studied for 
a number of polymer concentrations and for a 
polymer-colloid size ratio e=\ + 2Rg/d~ 1.8. 
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< [Piri] - Pstcp(zi)] [p{r2) - Pstop(^^2)] > , 

where p{r) is the colloidal density, rn — (x, y) is 
the direction parallel to the surface, and where 
we have deflned Pstc-p{z) = peQ{—z) + py 9(z) 
with 0{z) the Heavisidc function and pi^y the 
bulk density in the liquid and vapor region, re- 
spectively, where by "liquid" we mean the phase 
relatively rich in colloids and by "vapor" the 
phase relatively poor in colloids. Its Fourier 
transform is termed the surface structure factor 
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In Figure [H MC simulation results [J] for S{q) 
are shown for various values of the integration 
limit L. The Figure shows that the contribution 
to S{q) from short wavelength fluctuations (high 
q) increases with L. 

To analyze S{q), one needs to model the den- 
sity fluctuations in the interfacial region. In 
the capillary wave model (CW) [8], the fluctu- 
ating interface is described in terms of a two- 
dimensional surface height function /i(r||) 



The quantity studied in the simulations is the 
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FIG. 1: MC simulation results for the surface struc- 
ture factor (in units of d*) versus q (in units of 1/d) 
for various values of the integration limit L/W= 1, 
2, 3, 4 [4|. The dashed line is the capillary wave 
model. In this example e=1.8, r;p = 1.0. 

where Pq{z) — < p(f) >. In the extended capil- 
lary wave model (ECW), the exgansion in gra- 
dients of h{fn) is continued 



i3!,lM 



15|: 



pW = Po{z) - p'oiz) hifii) - 



Piiz) 



Ah{f\\) + . 



(4) 



The function pi{z) is identified as the correction 
to the density profile due to the curvature of the 
interface, A/i(r||) w— l/_Ri — I/R2, with i?i and 
i?2 the (principal) radii of curvature. 

With Eq.lg]) inserted into Eq.iH]), we find that 
S{q) equals the height-height correlation func- 
tion, S{q) — ShhiQ)i where 



5„^(g) = ydf|| e-'«-'-|i </i(fi,||)/i(f2,||)> . (5) 

Here we have assumed that the location of the 
interface, as described by the height function 
/i(n|), is given by the Gihhs equimolar surface 
[16| . which gives for po{z) and pi{z): 

dz[po{z) - pstep{z)] = 0, dz pi{z) =0. (6) 

Naturally, other choices are possible [5| and 
equally legitimate as long as they lead to a lo- 
cation of the dividing surface that is 'sensibly 

n . . 

coincident' [16] with the interfacial region. 



The height-height correlation function Shhiq) 
is determined by considering the free energy Afi 
associated with a surface fluctuation [8|, l9| . The 
inclusion of a curvature correction to the free 
energy is described by the Helfrich free energy 
[l7|. It gives for AO 

^^=\ Jj§^ ^(q) q' h{q) h{-q) , (7) 



with 



cr{q) = (T + kq^ 



(8) 



The coefficient k is identified as Helfrich's bend- 
ing rigidity [ITI . Il8| . It is important to realize 
that the bending rigidity, defined by Eqs.([7|) and 
([S]), depends on the choice made for the location 
of the dividing surface (here: the Gibbs equimo- 
lar surface for the colloid component). 

Using Eq.([7]), the height-height correlation 
function can be calculated |18j 
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(T{q) q^ a q^ + kq'^ + . . . 

Without bending rigidity (A: = 0) this is the clas- 
sical capillary wave result in the absence of grav- 
ity (dashed line in Figure [T]). When L is suf- 
ficiently large, the capillary wave model accu- 
rately describes the behavior of S{q) at low q. 

To model S{q) in the whole (/-range, we also 
include bulk-like fluctuations to the density: 



Pi.'r)^PQ(z)-p'o{z)h{r\\ 



Pi{z) 



A/i(f||)+<5pb(r). 
(10) 
Inserting Eq. (fTO| into Eq. ([1]) , one now finds that 



S{q) = Shh{q)+AfLSb{q). 



(11) 



The second term is derived from an integration 
into the bulk regions (to a distance L) of the 
bulk structure factor Sb{q) 

Sb{q) = l + Pb fdri2 e-'^l-^^'- [gir) - 1 ] . (12) 
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FIG. 2: MC simulation results J4| (circles) for the 
surface structure factor (in units of d*) versus q (in 
units of 1/d). The dotted line is the capillary wave 
model; the dashed line is the combination of the cap- 
illary wave model and the bulk correlation function; 
the solid line is the combination of the extended cap- 
illary wave model and the bulk correlation function. 
In this example e = 1.8, ??p= 1.0, L/W = 3. 

The density correlation function g{r) differs in 
either phase, but here we take for it gi{r) of the 
bulk liquid. This approximation may be justi- 
fied by arguing that close to the critical point 
there is no distinction between the two bulk cor- 
relation functions, whereas far from the critical 
point the contribution from the bulk vapor can 
be neglected since p„ « 0. The error is further 
reduced by fitting the L-dependent prefactor A/l 
to the limiting behavior of S{q) at qd ^oo. 

In Figured we show the result from Figure [1] 
for L/W— 3. For qd-^ 1 the results asymptoti- 
cally approach the capillary wave model (dotted 
line) . The dashed line is the result of adding the 
bulk-like fluctuations to the capillary waves: 
kBT 
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Figure [3 shows that Eq. p^ already matches the 
simulation results quite accurately except at in- 
termediate values oi q, qd~l. 

Finally, we include a bending rigidity in S{q): 
kBT 
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TABLE I: MC simulation results 4] for the polymer 
volume fraction rjp, liquid and vapor colloidal vol- 
ume fractions, -qi and i]^, surface tension a (in units 
of k-sT /d?), bending rigidity k (in units of fceT; in 
parenthesis the estimated error in the last digit), 
and yj —k/a (in units of d). 

The value of the bending rigidity is extracted 
from the behavior of S{q) at low q. The fact that 
the simulation results in Figure [5] are systemat- 
ically above the capillary wave model in this re- 
gion indicates that the bending rigidity thus ob- 
tained is negative, fc<0. Unfortunately, a nega- 
tive bending rigidity prohibits the use of S{q) in 
Eq. (|14p to fit the simulation results in the entire 
(/-range, since the denominator becomes zero at 
a certain value of q. It is therefore convenient to 
rewrite the expansion in q^ in Eq. (|14p as: 
knT , k 
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which is equivalent to Eq. (|14p to the order in 
q^ considered, but which has the advantage of 
being well-behaved in the entire g-range. The 
above form for S{q), with the bending rigidity 
used as an adjustable parameter, is plotted in 
Figure [5] as the solid line. Exceptionally good 
agreement with the MC simulations is now ob- 
tained for all q. In Table 1, we list the fitted 
values for the bending rigidity for a number of 
different polymer concentrations. 

Next, we investigate whether the value and 
behavior of k can be understood from a molec- 
ular theory. One should then consider a micro- 
scopic model for the free energy fJ to determine 
the density profiles po{z) and pi{z). Here, we 



consider the free energy density functional based 
on a squared- gradient expansion |7l. Il3l. Il4lll9j: 



n[p] = dr 



™|Vp(r-)p-|(Ap(f))2+g(p) 



(16) 



where the coefficients m and B are defined as 
m = -— jdri2 r'^U(r) , 

(dri2r^U{r). (17) 
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The integration over ri2 is restricted to the at- 
tractive part (r > d) of the interaction poten- 
tial U{r), for which we consider the Asakura- 
Oosawa-Vrij depletion interaction potential [10| : 



U{r) = 



-k^T 
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where the intermolecular distance is in the range 
1 < r/d < e. For explicit calculations, g{p) is 
taken to be of the Carnahan-Starling form: 



g{p) == UbT p\n{p)+kBT p 



(4r/ - 3772) 



-pp~ap , 
(19) 

where ry= (vr/G) pd^, /i = /^cocxj a-nd the van der 
Waals parameter a is given by 



dri2 U{r) . 



(20) 



The surface tension, to leading order in the 
squared-gradient expansion, can be determined 



from the usual expression 



A 

a 



cr = 2 ^Jm dp \/g{p) +p . 



(21) 



In the inset of Figure [3l the surface tension 
is shown as a function of the colloidal volume 
fraction difference, Ar] = ri£ — rjy. The squared- 
gradient ex pre ssion (solid line) is in satisfactory 
agreement |2(| with the MC simulations. 

The (planar) density profile pa{z) is deter- 
mined from minimizing the free energy func- 
tional Q[p\ in Eq. (|16p in planar symmetry. To 



also determine the density profile pi{z) from a 
minimization procedure, one should consider the 
energetically most favorable density profile for a 
given curvature of the surface. To set the cur- 
vature to a specific value, one adds to the free 
energy in Eg. lfTH)) an external field Vext{r) that 
acts a Lagrange multiplier. Different choices for 
Vextif) can then be made, but we choose it such 
that it acts only in the interfacial region: 



Vext{r) = Xp'aiz)Ahir« 



(22) 



with the Lagrange multiplier A set by the im- 
posed curvature. This choice for Vext{r) consti- 
tutes our fundamental 'Ansatz' for the determi- 
nation of pi(z). It improves on earlier choices 



made 



m 



2l| in the sense that the bulk den- 



sities are equal to those at coexistence and the 
density profile remains a continuous function. 

The minimization of the free energy, with the 
above external field added, using the fluctuat- 
ing density in Eq.Q yields the following Euler- 
Lagrange (EL) equations for Pq{z) and pi{z): 

g'ipo) = 2mp^'(z), 
g"ipo)piiz) - 2mp';{z) + Amp'a{z) (23) 
+2Bp'^'{z) + 2Xp',{z). 

The change in free energy Afi due a certain den- 
sity fluctuation is determined by inserting p(r) 
in Eq.(|4]) into the expression for il in Eq. lfTH)) . 
One finds that Afl is then given by the expres- 
sion in Eq.([71), with the bending rigidity [14 1 



k= -2m dz pi{z) Pq{z) - — jdz p'^{zf , 



(24) 



where we have used the EL equations in Eq. ([23|l . 

To determine po{z) we assume proximity to 

the critical point where g{p) takes on the usual 

double- well form. The solution of the Euler- 
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FIG. 3: Bending rigidity in units of ksT versus the 
volume fraction difference At;. The inset shows the 
surface tension in units of k^T/S. The solid lines 
are the gradient expansion approximation; filled cir- 
cles are the results from the MC simulations; the 
dashed line is the fit ^ —kja ~ 0.47 d. 



Lagrange equation in Eq. ((23|) then gives 



Q: 



Po(^) = \^P(^ Pv) - -J- tanh(z/2e) , (25) 

where ^ is a measure of the interfacial thick- 
ness which we shall define as i^ = r7i (Ap)^/(3 a), 
with the value of a given by Eq.([2T|l. To deter- 
mine pi{z) the differential equation in Eq. (j23|l is 
solved using the tanh- profile for po{z), yielding: 

3B Ap [l~ln(2cosh(z/2g))] 
IOto ^ cosh (z/2^) 

where we have used that A= — 2m + B/{5^'^). 

Inserting Eq. ((26|) into Eq. ((24| , one finds for k 

B{Apf _ Ba 



k = -- 



(27) 
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This expression indicates that the bending rigid- 
ity vanishes near the critical point with the same 
exponent as the surface tension, i.e. 

Ba 



k DC 



ex crd . 



(28) 



This scaling behavior should be contrasted to 
the usual assumption that fccxcr^^, i.e. that k 
approaches a finite, non-zero limit at the critical 

point \m,\2M- 

In Figure [3l the gradient expansion result in 
Eq. ([27|) for the bending rigidity is shown as the 



solid line. The bending rigidity is negative, in 
line with the simulation results, although the 
magnitude is significantly lower. 

To summarize, we have shown that to account 
for the simulated scattering function over the 
whole range of scattering vector q, including the 
intermediate range between low q (classical cap- 
illary wave theory) and high q (bulk-like fluctu- 
ations), one must take account of the interface's 
bending rigidity. Two of the important results 
are that the bending rigidity k for the interface 
between phase-separated colloid-polymer mix- 
tures is negative, and that on approach to the 
critical point it vanishes proportionally to the 
interfacial tension rather than, as had often been 
supposed, varying proportionally to the product 
of the tension and the square of the correlation 
length, thereby approaching a finite, non-zero 
limit. Both features of k are in accord with 
what is found in the simulations. The magni- 
tude of k obtained from the molecular theory is 



lower {y^^k/a ~ 0.13 d) than in the simulations 



[yj-k/ap^ 0.47 d; dashed line in FigurelS]). 
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